Convergence Bounds for Preconditioned Gmres Using Element-by-element Estimates of the Field of Values
نویسنده
چکیده
By combining element-by-element estimates for the field of values of a preconditioned matrix with GMRES-convergence estimates it is possible to derive an easily computable upper bound on the GMRES-residual norm. This method can be applied to general finite element systems, but the preconditioner has to be Hermitian and positive definite. The resulting upper bound for the GMRES-residual norm can be used to analyse a given preconditioner, or to optimize a parameter dependent preconditioner. In this paper we will apply this approach to derive a suitable shift for a so-called shifted Laplace preconditioner for the damped Helmholtz equation. Numerical experiments show that the shift that is derived in this way is close to optimal.
منابع مشابه
Field-of-Values Convergence Analysis of Augmented Lagrangian Preconditioners for the Linearized Navier-Stokes Problem
We study a block triangular preconditioner for finite element approximations of the linearized Navier–Stokes equations. The preconditioner is based on the augmented Lagrangian formulation of the problem and was introduced by the authors in [SIAM J. Sci. Comput., 28 (2006), pp. 2095–2113]. In this paper we prove field-of-values type estimates for the preconditioned system which lead to optimal c...
متن کاملPreconditioned Generalized Minimal Residual Method for Solving Fractional Advection-Diffusion Equation
Introduction Fractional differential equations (FDEs) have attracted much attention and have been widely used in the fields of finance, physics, image processing, and biology, etc. It is not always possible to find an analytical solution for such equations. The approximate solution or numerical scheme may be a good approach, particularly, the schemes in numerical linear algebra for solving ...
متن کاملField-of-values Analysis of Preconditioned Iterative Methods for Nonsymmetric Elliptic Problems
The convergence rate of Krylov subspace methods for the solution of nonsymmetric systems of linear equations, such as GMRES or FOM, is studied. Bounds on the convergence rate are presented which are based on the smallest real part of the eld of values of the coeecient matrix and of its inverse. Estimates for these quantities are available during the iteration from the underlying Arnoldi process...
متن کاملDomain decomposition preconditioning for high-frequency Helmholtz problems with absorption
In this paper we give new results on domain decomposition preconditioners for GMRES when computing piecewise-linear finite-element approximations of the Helmholtz equation −∆u − (k2 + iε)u = f , with absorption parameter ε ∈ R. Multigrid approximations of this equation with ε 6= 0 are commonly used as preconditioners for the pure Helmholtz case (ε = 0). However a rigorous theory for such (so-ca...
متن کاملSome Observations on the L2 Convergence of the Additive Schwartz Preconditioned GMRES Method ; CU-CS-865-98
Additive Schwarz preconditioned GMRES is a powerful method for solving large sparse linear systems of equations on parallel computers. The algorithm is often implemented in the Euclidean norm, or the discrete l norm, however, the optimal convergence theory is available only in the energy norm (or the equivalent Sobolev H norm). Very little progress has been made in the theoretical understanding...
متن کامل